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Abstract 

In the framework of the Lindblad theory for open quantum systems, a master equa- 
tion for the quantum harmonic oscillator interacting with a dissipative environment, in 
particular with a thermal bath, is derived for the case when the interaction is based 
on deformed algebra. The equations of motion for observables strongly depend on the 
iy"~j \ deformation function. The expectation values of the number operator and squared num- 

ber operator are calculated in the limit of a small deformation parameter for the case 
• of zero temperature of the thermal bath. The steady state solution of the equation for 

the density matrix in the number representation is obtained and its independence of the 
deformation is shown. 
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1 Introduction 

> ' 

• i-H , 

X ■ For more than a decade an interest has been consisting for the study of deformations of 
Lie algebras - so-called quantum algebras or quantum groups, whose structure produced 
important results and consequences in several branches of physics [TJ [2] . Their use in physics 
became intense with the introduction in 1989, by Biedenharn [3] and MacFarlane pE], of the 
g-deformed Heisenberg-Weyl algebra (g-deformed quantum harmonic oscillator). There are, 
at least, two properties which make g-deformed oscillators interesting objects for physics. The 
first is the fact that they naturally appear as the basic building blocks of completely integrable 
theories [5]. The second concerns the connection between g-deformation and nonlinearity. In 
Refs. [TJ [8] there have been considered some aspects of the physical nature of the q- 
oscillator and it was shown that it is subject to nonlinear vibrations with a dependence of 
the frequency on the amplitude. The g-deformed Bose distribution was obtained and shown 
how the g- nonlinearity produces a changed Planck distribution formula [6l EJ [9] . 



In the present paper we study the connection of quantum deformation and quantum 
dissipation, by setting a model for a harmonic oscillator interacting with a dissipative en- 
vironment, where the interaction between both is /- and g-deformed. The master equation 
which describes the dynamics of the oscillator represents a /-, respectively g-deformed version 
of the master equation obtained in the framework of the Lindblad theory for open quantum 
systems [10J. When the deformation vanishes, we recover the Lindblad master equation for 
the damped harmonic oscillator [LTJ[T2]. For a certain choice of the environment coefficients, 
a master equation for the damped deformed oscillator has also been derived by Mancini |13j . 
In Ref . [H] , Ellinas used the g-deformed oscillator for treatments of dissipation of a two-level 
atom and of a laser mode. In a previous paper [15] . using a variant of the Mancini's model 
[T3] , we derived a master equation for the /-deformed oscillator in the presence of a dissipative 
environment. 

The paper is organized as follows. In Sec. 2 we remind the basics about the /-deformed 
quantum oscillator, in particular the g-oscillator. In Sec. 3, using the Lindblad theory for 
open quantum systems, we derive a master equation for the harmonic oscillator in the presence 
of a dissipative environment, for the case when the interaction is deformed. The equations 
of motion obtained for different observables present a strong dependence on the deformation. 
In Sec. 4 we solve a linear system of two coupled first order differential equations for the 
expectation values of the number operator and squared number operator, in the limit of 
small deformation parameter and for zero temperature of the thermal bath. Then in Sec. 
5 we write in the number representation the equation for the density matrix and find the 
stationary state. In the particular case when the environment is a thermal bath we obtain 
the steady-state Boltzmann distribution of the ordinary harmonic oscillator, in which the 
deformation does not play any role. A summary is given in Sec. 6. 

2 /- and deformed quantum oscillators 

It is known that the ordinary operators {1, a, a', N} form the Lie algebra of the Heisenberg- 
Weyl group and the linear harmonic oscillator can be connected with the generators of the 
Heisenberg-Weyl Lie group. The /-deformed quantum oscillators [16J are defined as the 
algebra generated by the operators {1, A, A', N}. The /-deformed oscillator operators are 
given as follows [16] : 

A = af(N) = f(N + l)a, = /(JV)o+ = Jf(N + 1), (1) 

where iV = a^a. They satisfy the commutation relations 

[A,N]=A, [A\N] = -Ai (2) 
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and 

[A, A*] = (N + l)f(N + 1) - Nf 2 (N). (3) 

The function /, which is characteristics for the deformation, has a dependence on the defor- 
mation parameter such that when the deformation disappears, then / — > 1 and the usual 
algebra is recovered. Without loss of generality, / can be chosen real and nonnegative and 
it is reasonable from the physical point of view to assume [8] that /(0) = 1 and f(N) = 1 
for a suitable large N. The transformation (JTJ) from the operators a, a) to A, A* represents a 
nonlinear non-canonical transformation, since it does not preserve the commutation relation. 
On the other side, the transformation ([1]) does not change the dynamics, i. e. the equations 
of motion and, therefore, the transformation (JTJ) contains a symmetry for these equations. 
The notion of /-oscillators generalizes the notion of g-oscillators. Indeed, for 



where q is the deformation parameter (a dimensionless c-number), the operators A and A^ 
in Eqs. (JTJ) become the g-deformed boson annihilation and creation operators [3j H] . For the 
g-deformed harmonic oscillator the commutation relation ([3]) becomes 

[A,Ai] = [N + l}-[N]. (5) 

If q is real, then 

sinh(lng) 

and the condition of Hermitian conjugation (A*y = A is satisfied. In the limit q — > 1, q- 
deformed operators tend to the ordinary operators, lim 9 _ > i[A/'] = N and Eq. (J5J) goes to the 
usual boson commutation relation [A, A'] = 1. The g-deformed boson operators A and A* can 
be expressed in terms of the usual boson operators a and (satisfying [a,a^] = 1, N = a^a 
and [a, N] = a, [aJ, N] = —a') through the relations [TTJ [T8] : 



Using a nonlinear map [191, [20], the g-oscillator has been interpreted [6], [7] as a nonlinear 
oscillator with a special type of nonlinearity which classically corresponds to an energy de- 
pendence of the oscillator frequency. Other nonlinearities can also be introduced by making 
the frequency to depend, through the deformation function /, on other constants of motion, 
different from energy j6j [16]. Other examples of deformation functions / can be found in 
Refs. EH 12211231 EH. 
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3 Lindblad master equation with deformed dissipation 



Dissipative processes imply irreversibility, which means a preferred direction in time. There- 
fore, it is generally thought that a rigorous introduction of dissipation into a quantum system 
can be done in the framework of quantum dynamical semigroups [TQj, [25j [26] • The most gen- 
eral form of the generators of such semigroups was given by Lindblad |10j . This formalism 
has thoroughly been studied for the damped harmonic oscillator [HI [T2J, [27J, 12H] and applied 
to various physical phenomena, for instance, damping of collective modes in deep inelastic 
collisions in nuclear physics [12l [29] . 

According to the axiomatic theory of Lindblad [10J, the most general quantum Markovian 
master equation for the density operator p has the following form: 



l -[H,p{t)\ + ^Ea^PW,^] + IV^pW}}), (8) 



dp(t) 

dt " 1 2h 

j 

where H is the Hamiltonian of the system and the operators Vj and Vj , defined on the Hilbert 
space of H describe the interaction between system and environment. 

In this paper we consider that the open system is a harmonic oscillator with the Hamil- 
tonian 

1 2 mUj2 2 /n\ 

and the operators Vj are taken first degree polynomials in the basic observables p and q of 
the one-dimensional quantum system. Because in the linear space of these polynomials the 
operators p and q give a basis, there exist only two C-linear independent operators V\, V 2 
which have the form 

Vj = a 3 p + b jq , j = 1, 2, (10) 

with cij,bj complex numbers [27]. With these choices and introducing the annihilation and 
creation operators via the relations 



2muj 

the master equation takes the form 



q = J— — (a' + a), p = i\ /-^(a 1 - a), (11) 



-p = -~[H,p] + -{DAa ] a ] p + pa)a ] - 2a 1 pat) 
dt n 2 

+(D 2 + A)(apa t - a ] ap) + (D 2 - A)(a t pa - paa ] ) + H.c.}, (12) 
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where H = {l/2)huj{aa [ + a) a) and 

Di = \(mu)D qq - ^ + 2iD pq ), D 2 = \{muD qq + ^). (13) 
n muj n moo 

We used the notations: 

D m = 9 J2 Kf> D pp = 9 J2 l 6 i! 2 ' D pi = D ip = -| Re E a *j b ii X = ~ lm a ] h jl U ) 

Z j=l,2 Z j=l,2 Z j=l,2 j=l,2 

where D pp , D qq and D pq are the diffusion coefficients and A the dissipation constant. These 
quantities satisfy the following fundamental constraints [Tlj [12] : 

h 2 X 2 

i)D pp >0, n)D qq >0, iii) D pp D qq - D 2 pq > —. (15) 
In the particular case when the asymptotic state is a Gibbs state 



H , H 



Pg(oo) = e ^/Tre (16) 

the diffusion coefficients reduce to 

A hu A h hu 

D pp = -hrru, coth—, D qq = -—cot\i—, D pq = 0, (17) 

where T is the temperature of the thermal bath and k the Boltzmann constant. The relations 
(1T5]) ensure the positivity of the density operator. Due to the dissipative terms in Eq. (TT2"|) . 
the energy of the open system is dissipated into the environment, but noise arises also (in 
particular, thermal noise), since the environment also distributes some of its energy back to 
the system. 

Next we propose a model where we deform the dissipative part in Eq. ( TT2l) by replacing 
the operators a and by the /-deformed operators ([1]). Then the master equation ([121 
becomes: 

+ l -D 1 {}\N - l)f(N)aWp + paWf(N + l)f(N + 2) - 2f(N)a^pa^f(N + 1)} 

+ l -D\{f{N + l)f(N + 2)aap + paaf(N - l)f(N) - 2f(N + l)apaf(N)} 

- 1 -{D 2 + \){Nf 2 (N)p + pNf 2 (N) - 2f(N + l)apa^f(N + 1)} 

- 1 -{D 2 -\){(N+ l)f 2 (N + l)p + p(N + l)f 2 (N + 1) - 2/(JV)oW(JV)}. (18) 

We stress that in the proposed model we do not introduce the deformation into the com- 
mutator containing the oscillator Hamiltonian H, but we do it only in the dissipative part 
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of the master equation, which describes the influence of the environment on the oscillator. 
The deformation of the oscillator Hamiltonian would imply the deformation of the diffusion 
coefficients too and this model has been treated in our previous paper [15J. Obviously, the 
obtained master equation preserves the Hermiticity property of the density operator and the 
normalization (unit trace) at all times (if at the initial time it has these properties). In the 
limit / — > 1 the deformation disappears and Eq. ( fl8l) becomes the Markovian master equation 
for the usual damped harmonic oscillator, obtained in the Lindblad theory for open quantum 
systems. In the particular case of a thermal equilibrium of the bath at temperature T with 
diffusion coefficients ( 1T7|) . the master equation ( flBl) takes the form 

-^A(coth ^ + l){Nf(N)p + pNf(N) - 2f(N + l)apctf(N + 1)} 
-7^(coth ^ - 1){(N + l)f(N + l)p + p(N + l)f(N + 1) - 2f(N)tJpaf(N)}. (19) 
If the bath temperature is T = 0, the master equation ffl9|) simplifies and becomes 

^ = -iu[N, p] - \{Nf 2 (N)p + pNf 2 (N) - 2f(N + l)ap<J f(N + 1)}. (20) 

4 Equations of motion for < N > and < TV 2 > 

The meaning of the master equation becomes clear when we transform it into equations 
satisfied by various expectation values of observables involved in the master equation, < 
B >= Tr[p(t)B], where B is the operator corresponding to such an observable. We give an 
example, multiplying both sides of Eq. ( fl9l by the number operator N and taking the trace. 
Then we obtain the following equation of motion for the expectation value of iV: 

j t < N >= A{(coth ^ - 1) < (N + l)f\N + 1) > -(coth ^ + 1) < Nf(N) >}. (21) 

We consider another example, namely for the operator N 2 and we obtain the following equa- 
tion: 

j t < N 2 >= A{(coth ^ - 1) < (2JV + l)(N + l)f 2 (N + 1) > 

-(coth ^ + 1) < (2N - l)Nf(N) >}. (22) 

We observe that the equations of motion for the expectation values contain complicated 
nonlinearities introduced by the deformation function / and, in addition, they do not form a 
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closed system. In order to introduce some procedure for truncating the set of these coupled 
equations, we consider the simplest case of deformation - the g-deformation. In this case, 
using Eq. ([4]), Eqs. ( |2TI) and f l22l) take the following form: 

j t < N >= A{(coth J£L-i) < [N+l]> -(coth + 1) < [N] >}, (23) 

I < N 2 >= A{(coth ^ - 1) < (2N + 1)[N + 1] > 

-(coth^ + l)<(2iV-l)[iV]>}. (24) 

Now we write the g-operator [N] (JS]) in the lowest order of the Taylor expansion in the small 
deformation parameter r (we consider the case of real q = exp r) [2] : 

\N} = N (iV-iV 3 ). (25) 

6 

With this expression, Eqs. ff23|) and ff24|) become: 

^- < N >= Aicoth -44 — 1 — 2 < iV > 
dt 1 2fcT 

+ — [3 coth—- < MiV + 1) > - <2A^ 3 + 3iV 2 + A^>]}, (26) 
6 2kT 

— < N 2 >= Aicoth 44 < 4iV + 1 > - < AN 2 + 2iV + 1 > 
at 2kT 

+4f coth t44 < 8A^ 3 + 9iV 2 + > - < 4N 4 + 6N 3 + 5N 2 + 3N >]}. (27) 
6 2fcT 

Now we make use of the relation 

N 3 = a t3 a 3 + 3N 2 - 2N (28) 

and make the assumption that < a^ 3 a 3 > can be neglected for not too high values of the 
number operator N. Then, replacing A^ 3 by 3A^ 2 — 2N in Eqs. (1261) and (1271) . we obtain: 

< jV >= A {—coth 3) < AT 2 > 

dt 1 2 K 2kT ' 

t 2 fiuj , fou; . . 

_[2 - - (coth — + 1)] < iV > + coth — - 1}, (29) 

4 < N 2 >= A{[— (11 coth 44 " 17) - 4] < N 2 > 
dt U 2 V 2kT ' J 

, huj t 2 . Tilo , fru; . 

+ 4 coth— — -2 5 coth— — -11 < JV > + coth — — - 1}. 30 

L 2kT 2 v 2kT n 2kT s v ; 



By this truncation, we obtained a closed non-homogeneous linear system of first order coupled 
differential equations for the expectation values < N > and < iV 2 >, which can analytically 
be solved. For simplicity, we will solve this system in the case of the temperature T = of 
the thermal bath. Then Eqs. (1291) and (|30|) become: 

— < N >= -r 2 A < iV 2 > +(r 2 - 2)A < iV >, (31) 
< N 2 >= -(3r 2 + 4)A <N 2 > +(2 + 3r 2 )A <N> . (32) 



The integration of this system is straightforward. Introducing the vector 

< N(t) > 

< N 2 (t) > 



and the matrix 

M = \ l 

,2 + 3r 2 -3r 2 -4 



M = \[f lQ \ J' J, (3D 



the system ( 13~TT) - ( 1321) can be written 

dS(t) 



MS(t). (35) 



dt 

Now M can be written as M = R~ l FR with F the diagonal matrix that has as elements the 
roots of the characteristic equation of the linear system. It follows that the solution of the 
system is 

S{t) = R- l e Ft RS{0). (36) 

We obtain 

p-i Ft p - 1 / (2 + 3r 2 )e- 2A * - r 2 e - 2A ( 2 +- 2 )* -r 2 (e" 2A ' - e^ 2+T ^) \ . . 
K e K ~ 2 (l + r 2 ) v(2 + 3r 2 )(e- 2A * - e - 2A < 2+ ^) -r 2 e - 2A * + (2 + 3r 2 ) e - 2A ( 2+ ^V [6 

and the solutions are the following: 



< N(t) >= 1 {-r 2 (e- 2Xt - e - 2A ( 2 +- 2 )*) < N 2 (0) > 

2(1 + T ) 

+ [(2 + 3r 2 )e- 2A * - r V 2A ( 2 +- 2 )*] < iV(0) >}, (38) 



< N 2 (t) >= o/i ' - {[-r 2 e' 2Xt + (2 + 3r 2 ) e - 2A ( 2+T2 ) t ] < iV 2 (0) > 



1 

2(l + r 2 

+ (2 + 3r 2 )(e- 2A * - e - 2A(2+r2)t ) < N(0) >}, (39) 
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or, in a more simple form 

1 



< N{t) >= - 2 A e~ 2Xt [-T 2 < N 2 (0) > +(2 + 3r 2 ) < JV(0) >] 
2(1 + r ) 



2 e -2A(2+^) t(< jv(0) > _ < iy2(0) >)} ^ (4Q) 



-r 



< N 2 (t) >= 1 {[e~ 2Af [-r 2 < iV 2 (0) > +(2 + 3r 2 ) < JV(0) >] 
2(1 + r ) 

+ (2 + 3r 2 )e- 2A{2+r2)t (< iV 2 (0) > - < iV(0) >)}. (41) 
Keeping only the relevant terms in small r 2 , we can write the solutions as follows: 

< N(t) >= e~ iXt — (< iV 2 (0) > - < N(0) >) 

-e" 2A '[y < iV 2 (0) > -(1 + y) < iV( ) >], (42) 
< N 2 (t) >= e~ 4Xt (l + y - 2Ar 2 t)(< iV 2 (0) > - < iV(0) >) 

2 2 

" e ~ 2A< [y < ^ 2 (0) > "(1 + — ) < N(0) >). (43) 

In the case when q is a phase (q = exp(zr),r real), the expectation values of the number 
operator and squared number operator are obtained by making the replacement r 2 -> -t 2 
in the corresponding expressions (T42|) and (T43j) . In the limit of long times we obtain < 
N 2 (oo) >=< N(oo) >= 0, like in the case of the undeformed case. This is as we expected, 
since we are in the case of T = 0. In the limit of vanishing deformation (r — > 0) we obtain: 

< N(t) >= e' 2Xt < N(0) >, (44) 



< N 2 (t) >= e- 4Xt {< N 2 (0) >-< N(0) >) + e- 2Xt < N(0) >, (45) 

in concordance with the previous results obtained in the study of the damped harmonic 
oscillator in the Lindblad theory [30]. The exponential decreasing of the expectation values 
of the number operator and of the squared number operator are represented graphically in 
Fig. 1 for both the cases when q = expr and q = exp(zr), for a small real deformation 
parameter r = y/0.2. For comparison we represented also the evolution of these observables 
in the case when the deformation is absent. The equations of motion for other observables 
contain also nonlinearities introduced by the deformation function /, like in the considered 
case. 
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At 



0.2 0.4 0.6 0.8 1 



Figure 1: Dependence of (a) averaged number operator < N > given by Eq. (|42p and of (b) 
averaged squared number operator < N 2 > given by Eq. (j43p on the dimensionless time Xt, 
for < N(0) >= 3, < iV 2 (0) >= 9 and r 2 = 0.2. The cases of q — expr, q = exp(zr) with r 
real and q = are represented by solid, dashed and dotted lines, respectively. 
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5 Equation for the density matrix 



Let us rewrite the master equation (fT8l) for the density matrix by means of the number 
representation. Namely, we take the matrix elements of each term between different number 



states denoted by \n >, and using N\n >= n\n >, a + \n >= y/n + l\n + 1 > and a\n >- 
y/n\n — 1 >, we get 

— — = -iu(m - n)p mn 

- 1 -{{D 2 + \)(mf(m) + nf(n)) + (D 2 - A)((m + l)f(m + 1) + (n + l)f(n + l))}p mn 

+ (D 2 + X)y/(m + l)(n + l)f(m + l)f{n + l)p m+ i, n +i 
+(D 2 - A)v / rnn/(m)/(n)p m _i )n _i 
-Diy/m(n + l)f(m)f(n + l)p m -i, n +i 
-D{^Jm + l)nf(m + 1) / (n) p m +i,n-i 



+ -D^(n + l)(n + 2)f(n + l)/(n + 2)p m , n+2 

+ ^Di^m(m - l)/(m - l)/(m)p TO _ 2 ,„ 

+ ^D 1 * v /(m + l)(m + 2)/(m + l)/(m + 2)p m+2 , n 

+^ lV /n(n - l)/(n)p m , n _ 2 .(46) 

Here we used the abbreviated notation p mn =< m\p{t)\n > . This equation, complicated 
in form and in indices involved, comprises not only the density matrix elements in sym- 
metrical forms, such as p m ±i,n±i, but also those matrix elements in asymmetrical forms like 
Pm±2,n, Pm,n±2 and p m =Fi,n±i- Eq. (HBT) gives an infinite hierarchy of coupled equations for the 
off-diagonal matrix elements. Their evolution in time depends on the diffusion and friction 
coefficients and on the deformation. When Di = 0, the diagonal elements are coupled only 
amongst themselves and not coupled to the off-diagonal elements. In this case the diagonal 
elements (populations) satisfy a simpler equation: 

**M = _{(D 2 _ A)(n + l)f\n + 1) + (D 2 + \)nf{n)}P{n) 
at 

+ (D 2 + A)(n + l)f{n + l)P(n + 1) + {D 2 - \)nf{n)P(n - 1), (47) 
where we have set P{ri) = p nn . For a g-deformation, Eq. (I47p takes the form 

^ = -{(D 2 - A)[n + 1] + (D 2 + X)[n]}P(n) 
+(D 2 + X)[n + l]P(n + 1) + (D 2 - A) [n]P(n - 1). (48) 
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We define the transition probabilities 

t+(n) = (D 2 -\)(n + l)f(n + l), i_(n) = (D 2 + X)nf 2 (n), (49) 
which for a g-deformation look like 

t+ (n) = (£> 2 - A) [n + 1] , f _ (n) = (Z> 2 + A) [n] . (50) 
With these notations Eq. (j4"7j) becomes: 

= t + (n - l)P(n - 1) + t_(n + l)P{n + 1) - [t+(n) + t_(n)]P(n). (51) 
The steady state solution of Eq. (1511) is found to be 

= (52) 

We remark that in the steady state the detailed balance condition holds: 

t_(n)P{n) = t+{n - l)P(n - 1). (53) 

In the particular case of a thermal state, when the diffusion coefficients have the form (ITT)) , 
the stationary solution of Eq. (j5Tl) takes the following form: 

n) = Z- 1 exp[-^(2n + l)] ) (54) 



with 



Z-W^exp^, (55) 



Z being the partition function for the ordinary harmonic oscillator 



00 hi s 1 

n=0 2fcT 2 Smh 2fcT 

The distribution (1541) can also be written 

P^H = Z- 1 exp(-^), (57) 

where E n = hu(n + 1/2) is the ordinary expression for the energy of the harmonic oscillator. 
Expressions (1541) and (1571) have the known form of the Boltzmann distribution. We notice 
that in the expressions (152]) and (1541) of the stationary distribution the deformation of the 
dissipation does not play any role, as expected, due to the fact that in the considered model 
the deformation is introduced only in the interaction of the system with the environment and 
the Hamiltonian of the system is kept undeformed. 
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6 Summary 



The Lindblad theory of open systems based on quantum dynamical semigroups provides a 
selfconsistent treatment of damping as a possible extension of quantum mechanics to open 
systems. Our purpose was to study the dynamics of the quantum harmonic oscillator in 
interaction with a dissipative environment, with a /- and g-deformed type of interaction. We 
have deformed only the operators in the dissipative part of the Lindblad equation, but not 
the Hamiltonian itself. In this case the diffusion and dissipation coefficients which model the 
influence of the environment on the oscillator are constant quantities for a given temperature 
and do not depend on the introduced nonlinearities (number operator). In the limit / — > 1, 
when deformation disappears, the obtained master equation for the damped oscillator with 
deformed dissipation becomes the usual master equation for the damped oscillator obtained 
in the framework of the Lindblad theory. The equations of motion for the observables of 
the considered system are strongly nonlinear and they form a set of coupled differential 
equations which is not closed. By truncating the set of coupled differential equations for the 
expectation values of the number operator and squared number operator, we were able to 
obtain a closed system of coupled differential equations in the case of a small g-deformation. 
The obtained solutions show explicitly the dependence on the small deformation parameter 
t. We have also derived the equation for the density matrix in the number representation 
and in the case of a thermal bath we obtained the stationary solution, which is of the form 
of the ordinary Boltzmann distribution, as expected, since we did not deform the harmonic 
oscillator Hamiltonian. 

A more complicated model one obtains by introducing a /-deformation for the harmonic 
oscillator operators, both in the Hamiltonian of the system and in the dissipative part of the 
master equation. Work is in progress in this direction. 

The master equation for the damped harmonic oscillator with deformed dissipation is an 
operator equation and it could be useful to study its consequences by transforming it into more 
familiar forms, such as the partial differential equations of Fokker-Planck type for the Glauber, 
antinormal ordering and Wigner quasiprobability distributions or for analogous deformed 
quasiprobabilities [8] associated with the density operator. It could also be interesting to find 
the states which minimize the rate of entropy production, which play an important role in 
the description of the decoherence induced by the environment. 
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